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ABSTRACT 

The Kepler Mission seeks to detect Earth-size planets transiting solar-hke stars in its ~115 deg^ 
field of view over the course of its 3.5 year primary mission by monitoring the brightness of each 
of ~156,000 Long Cadence stellar targets with a time resolution of 29.4 minutes. We discuss the 
photometric precision achieved on timescales relevant to transit detection for data obtained in the 33.5- 
day long Quarter 1 (Ql) observations that ended 2009 June 15. The lower envelope of the photometric 
precision obtained at various timescales is consistent with expected random noise sources, indicating 
that Kepler has the capability to fulfill its mission. The Kepler light curves exhibit high precision 
over a large dynamic range, which will surely permit their use for a large variety of investigations in 
addition to finding and characterizing planets. We discuss the temporal characteristics of both the 
raw flux time series and the systematic error-corrected flux time series produced by the Kepler Science 
Pipeline, and give examples illustrating Kepler's large dynamic range and the variety of light curves 
obtained from the Ql observations. 

Subject headings: techniques: photometric — methods: data analysis 



1. INTRODUCTION 

The Kepler Mission is designed to be capable of find- 
ing ~50 Earth-like planets transiting solar-like stars 
in the solar neighbor hood if such planets are common 
(|Borucki et al. 1 120101 ) . Specifically, Kepler is required 
to obtain a signal-to-noise ratio (S/N) of 4 cr for an 
84 ppm deep, 6.5-hr transit of a G2V star wit h a Ke- 
pler magnitude {Kp) of 12 (|Koch et al. IIMot ). This 
implies that the rms noise on 6.5-hr intervals must be 
below 20 ppm, including contributions from stellar vari- 
ability. Since 2009 May 12, Kepler has been observing 
^156,000 stars at 29.4 minute intervals as Long Cadence 
(LC) targets for the primary purpose of detecting tran- 
siting planets. Of these targets, 512 are also sampled at 
1 minute intervals to suppo rt asteroseismic charact eriza- 
tion of planetary host stars (Gill iland et al. Il2010al lbl). as 
well as precise timing of transits. 

In Sectoin[2]we examine the photometric precision ob- 
tained to date on the first 33.5 days of data obtained 
during Quarter 1 (Ql) which began May 13 and ended 
June 150. We discuss temporal characteristics of the raw 
and systematic error-corrected flux time series in Section 
[3l and present sample light curves in Section|4l Summary 
remarks appear in Section [5l 

2. PHOTOMETRIC PRECISION 

|Jon .Jenkins@nas a.gov| 

Kepler's science operations are broken into quarterly seg- 
ments, each of which is ~93 days long and is subdivided into three 
monthly segments. Data are downlinked at these monthly intervals 
HHaas et al. l[20Toh . 



To address the question of the photometric precision 
obtained from the Kepler LC data, we examined system- 
atic error-corrected flux time series furnished by the Pre- 
search Data Conditioning (PDC) component of the Sci- 
ence Operations Center (SOC) pipelineO, hereafter called 
PDC time series. At the end of each quarter, the space- 
craft rolls 90° about the boresight in order to reorient 
the solar arrays toward the Sun, and the target stars ro- 
tate onto new det ectors in the Focal Plane Array (EPA; 
iHaas et al. Il2010f ). Each of these quarterly segments is 
processed by first calibrating the raw pixels and then 
the Photometric Analysis (PA) component subtracts the 
sky background and extracts simple aperture photomet- 
ric brightness estimates for each star. Systematic errors 
and instrumental effects are identified and removed from 
the raw flux time series by PDC. 

The PDC flux time series for all 156,097 stars were pre- 
pared for analysis by subtracting the median flux value, 
normalizing by the same, and then removing a fltted 4th 
order polynomial across all 1639 LC samples. The collec- 
tion of standard deviations of the normalized, detrended 
flux time series represents the photometric precision on 
29.4 minute intervals as a function of stellar magnitude, 
and this is plotted in panel A of Figure [TJ 

The standard propagation of errors is applied to each 
data processing step so that flux uncertainties are avail- 
able along with the flux measurements themselves. These 
uncertainties include shot noise from photon counting 
statistics from the direct stellar flux, sky background 

^ Details of the SOC science pipeline are discussed in 
IJenkins et"al~l (l20TOfl . 
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and shutterless readout smear, terms such as read noise 
and quantization noise, and the effects of processing the 
data, for example, pixel-level calibrations and sky back- 
ground removal. The uncertainties depend on the de- 
tails of the size of the photometric aperture, the Pixel 
Response Function, and other considerations. It is im- 
portant to note that these uncertainties do not include 
any contributions from intrinsic stellar variability. 

Plotting the individual uncertainty estimates 
would clutter the figure, so we fit a simple an- 
alytic expression "by eye" to the lower enve- 
lope of these uncertainties at the LC timescale: 
o'lowcr = y/c + 7 X 10^ max(l, Kp/14:)^/c where 
c = 3.46 X io[0-4*(i2-ifp)+8] is the number of de- 
tected per LC sample, and Kp is the Kepler Input 
Catalog (KIC) magnitude for each star. The component 
c is a Poisson term while the {Kp/IA)'^ accounts heuris- 
tically for the increased importance of sky background 
and other noise sources for the dimmer stars with 
Kp > 14. We also fit "by eye" an analytic expression 
to the upper envelope of this distribution as a function 
of Kp: (Tuppor = Vc -f 7 X lO'^/c. In this case a simple 
constant background term was sufficient. The lower and 
upper uncertainty envelopes appear in Figure [1] as solid 
and dashed green curves, respectively. 

Three populations of stars can be seen in panel A. We 
take stars with KIC entries of logg < 4 to be giants 
or sub-giantfl For the most part, the main sequence 
stars float just above and follow the profile of (Tiower as 
a dense "finger" of points extending over the full mag- 
nitude range of the stars. The fractional width of this 
"finger" is approximately 1.5 near Kp = 12 and flares 
out to a factor of 2 at Kp = 16. This is not unex- 
pected as the dim stars are more sensitive to background 
sky noise, point spread function variations, and detec- 
tor properties, which vary significantly from detector to 
detector. A high "cloud" of stars, composed of both gi- 
ants and dwarfs, hovers from 3,000 to 10,000 ppm across 
nearly the full magnitude range. These are variable stars 
of all flavors, and possibly some new ones observed for 
the first time. There is a sharp change in t he density 
of points for Kp > 14 (jBatalha et al. I l2010f ). This in- 
cludes the fact that no stars classified as giants with 
Kp > 14 were selected. Giant stars cluster on a hori- 
zontal band centered between ~400 and ~700 ppm with 
Kp < 14. The fact that giant stars exhibit significantly 
more photometric variability on timescales of minutes 
to hours co mpa red to dwa rf stars has been previously 
observed bv ICilliland I (|200S l in Hubble Space Telescope 
observations of bulge stars, and more recently, by CoRot 
(jAigrain et al. 1 120091 ). Giant stars universally appear to 
exhibit solar-like oscillations, and their increased size rel- 
ative to solar-type stars causes these oscillations to occur 
at much lower frequencies a nd larger amplit udes than 
observed for solar-type stars (|Gilliland et al. I i2010b). A 
total of 350 points lie below the aiower- These stars suf- 
fer from excess flux in their apertures likely from nearby 
brighter stars or bleeding charge from saturated stars in 
the same column. When these points are plotted at their 

^ The KIC is 90% reliable for the purposes of separating giants 
from dwarf stars for Kp < 13 a nd the corresponding uncertainties 
in \ogg are at the 0.5 dex level (|Batalha et al. II2010I V 



instrumental magnitudes (i.e., their measured counts are 
converted to magnitudes via the inverse of the relation 
given above), rather than by catalog magnitude, Kp^ 
they are all above 

As the central transit of an Earth analog lasts ^13 hr, 
we examine the photometric precision relevant for detect- 
ing Earth-like transits lasting 6 hr or more. Each PDC 
flux time series considered was first high pass filtered to 
remove variations on timescales of 2 days and longer and 
then convolved with a 6-hr wide "boxcar" averaging fil- 
ter. The high pass filtering was accomplished by fitting 
a 2nd order polynomial to a sliding analysis window 2 
days wide about each point and removing the resulting 
polynomial at the central point for each window loca- 
tion. We calculated the standard deviations of the time 
series resulting from this process as estimates of the 6 hr 
precision of these data and present these in panel B of 
Figure [TJ The LC uncertainty envelopes have been nor- 
malized by ^/l2 to account for square root statistics. The 
giant stars (red points) occupy a band between 200 and 
500 ppm. The main sequence stars continue to lie in a 
dense "finger" just above CTiowor/V^- The bottom enve- 
lope of the main sequence population tracks the expected 
measurement uncertainties quite well. A total of ~2000 
points lie below the o'lowcr/'v/lSj but since there are fewer 
independent measurements at the 6 hr level, this is ex- 
pected. 90% of the 2,665 stars with 11.5 < Kp < 12.5 
in the dense main sequence "finger" are below 30 ppm, 
65% are below 25 ppm, 50% are below 23 ppm, and 31% 
are below 21 ppm (which corresponds to 20 ppm at 6.5 
hours). 

Kepler is achieving extremely high precision over a re- 
markably large dynamic range considering that the mis- 
sion design was driven largely by the 20 ppm, 6.5-hr pre- 
cision requirement for Kp = 12 G2V stars. Clearly, stars 
do exhibit stellar variability, but there are stars that fall 
near the aiowcv/V^ down to at least Kp = 9 and there 
is a Kp = 8 star exhibiting noise on 6 hr intervals of 4 
ppm! We believe that further improvements in photo- 
metric precision will be made as we gain experience with 
the spacecraft and the photometric data. 

3. TEMPORAL CHARACTERISTICS OF LC LIGHT CURVES 

Although the photometric precision is a key consid- 
eration in assessing Kepler's ability to detect transiting 
terrestrial planets, it is not the only one. The precision 
metrics discussed in Section [2] only speak to the mean 
noise power per 29.4 minuts or per 6 hr interval. Here 
we discuss temporal characteristics of the raw and PDC 
light curves of importance for transit detection and for 
non-exoplanet investigations. 

There are systematic effects that occur on timescales 
much longer than those relevant for detecting transits. 
These include Differential Velocity Aberration (DVA), 
which causes stars to move as much as 0.6 pixels over a 93 
day period due to the large size of the field of view (FOV). 
Kepler has also experienced secular pointing drifts of ~ 
0.05 pixels over the Ql observations, and there have been 
thermal systematic effects occurring on ti mescales of sev- 
eral d ays and longer seasonal variations (jJenkins et al. I 
120101 ) as a consequence of the Sun's sky position rotating 
about the telescope by about 1° day~^. A bulk change 
in focus (including tip-tilt terms as well as a "piston" 
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term) occurred in part due to desorption of water from 
the telescope structure and seasonal thermal effects. The 
gradual change in focus has caused the measured flux of 
stars on the outer periphery of the FOV to drop by about 
0.1% over Ql while stars toward the interior (on the op- 
posite side of best focus) have experienced corresponding 
increases in flux. 

Systematic errors on these timescales are apparent es- 
pecially in the raw photometric data, but are largely 
suppressed in the PDC time series. Astrophysical sig- 
natures occurring on timescales similar to the long term 
seasonal variations in Kepler's systematic errors will in- 
evitably be attenuated by PDC. We do not intention- 
ally remove astrophysical signatures on any time scale. 
However, the conditioning of the data simply does not 
allow us to distinguish between astrophysical signatures 
and systematic errors if the former are correlated with 
ancillary engineering data and science data diagnostics 
used in the systematic error removal process. Neverthe- 
less, the PDC time series display a remarkable wealth 
of astrophysics at virtually all timescales probed by the 
observations. 

3.1. Pointing Drift and Focus Changes 

As reported in iJenkins et al. I ()2010l ) the photometer 
is quite sensitive even to the slightest change in the ther- 
mal environment. Figure |5] shows changes in a) the focus 
distance local to module 2.1, b) mean row, and c) mean 
column positions of the ensemble of stars located on mod- 
ule 2.1 near the outside corner of the FPA over a 15 day 
interval during Ql. These time series have been filtered 
to remove trends on timescales longer than 2 days, which 
are not relevant for transit detection. The focus distance 
of this module appears to be changing by about 1 /im ev- 
ery time a heater on one reaction wheel assembly turns 
on to keep the wheels above the setpoint temperature. 
These are the downward spikes in the focus distance, 
which represents the change in the relative distances be- 
tween the stars over the observations. The mean row and 
column positions of these stars also experience changes in 
response to this heater cycle. Additionally, they experi- 
enced excursions every 1.7 days lasting about 8 hr due to 
an eclipsing binary on the Fine Guidance Sensor (FGS) 
target list that was modulatin g the spacecraft pointing 
by '^l mpix during eclipses (H aas et al. Il2010l l. 

Similar features appear in the raw flux time series of 
many LC targets. Some respond to the pointing excur- 
sions, while others do not. Most respond to focus changes 
at some level. Fortunately, PDC was designed to remove 
systematic effects of this kind and does a good job of re- 
moving such signatures for relatively quiet stars and for 
stars that exhibit coherent oscillations with slowly vary- 
ing frequency components. The harmonic components of 
the latter type star are removed prior to removing sig- 
natures correlated with information such as presented in 
Figure [H and are then restored at the end of the cor- 
rection process. There are incoherently variable targets 
for which it is difficult to separate the intrinsic stellar 
variability from systematic effects. PDC detects when 
this is the case and provides a much simpler detrending 
scheme for such targets. Figure |3] exhibits one example 
of raw and PDC flux time series for one star located on 
module 2.1 which is one of the modules most sensitive to 
the short-term thermal effects. 



We are making excellent progress with respect to un- 
derstanding the behavior of the photometer. However, 
it is certain that much future work will consist of mon- 
itoring the systematics available through ancillary engi- 
neering data and science data diagnostics so that we can 
best remove the effects from the flux time series while 
retaining as much intrinsic stellar variability as possible. 

3.2. Cosmic Rays 

Pre-launch expectations for the cosmic rays Kepler 
would encounter on orbit were informed by a radiation 
model developed by Ball Aerospace Technologies Cor- 
poration (Neal Nickles, private communication 2003). 
Based on this model, we adopted a conservative hit rate 
of 5 cm~^ sec~^, and the mode of the distribution of to- 
tal charge deposited per hit was expected to be 2500 e^. 
These facts translate into an expectation that every pixel 
will receive a direct cosmic-ray hit rate of ~ 3 day"^ and 
that each such event would contribute charge to ~4 ad- 
jacent pixels. Analysis of dark frames acquired prior to 
releasing Kepler's dust co ver indicate a cosmic ray depo - 
sition rate of 10 cm'^ s'^ ()Van Cleve fc Caldwell Il2009l ) . 
consistent with our pre-launch expectations. Thus, cos- 
mic rays do not pose a significant noise problem for tran- 
sit detection. 

3.3. Image Artifacts 

Image artifacts introduced by some circuits in the ana- 
log electronics for Kepler's Focal Plane Array were iden - 
tified during pre- flight testing (jCaldwefl et ah 1 120101 ). 
Although a small fraction of the FOV is impacted at any 
one time, the effects are either amenable to correction in 
data processing or can be identified and excluded from 
consideration at little risk to the primary mission goals. 
The risk to the hardware and cost of attempting a fix 
to the electronics were too high to proceed with a direct 
intervention with the flight hardware. The Kepler Sci- 
ence Office and SOC are engaged in developing, testing 
and implemen ting software m itigations for these image 
artifacts (Jen kins et al. I [20101 ). Typical amplitudes of 
these artifacts are < 1 ADU pixel" ^ read~^, so that in 
most cases their presence in the flux time series of stars 
can be detected only by sophisticated analyses. Once 
the mitigations are in place we expect to recover the re- 
quired photometric pe rformance for over 90% of the FOV 
(ICaldwell et al. imuh . 

3.4. CCD Defects 

A few otherwise good targets land on column or pixel 
defects. Secular drift in pointing together with DVA car- 
ries stars by up to ^0.7 pixels over each quarter, and 
as an affected star moves across a column defect, short- 
term pointing errors can lead to abrupt changes in the 
measured flux value so that the flux appears to "toggle" 
between two states for a period of time. Currently we 
make no attempt to identify and correct this behavior. 
Panels A and B of Figure [3] show the light curve for KID 
5286794 illustrating this behavior. Very few targets are 
affected by this phenome non, given the paucity of CCD 
defects, as documented in ICaldwell et al. 1 ()2010D . 

4. EXAMPLES OF LIGHT CURVES 

In this section we present a potpourri of light curves 
obtained during Ql observations that illustrate Kepler's 
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stunning photometric precision and the remarkable di- 
versity of stehar photometric behavior across a wide 
range of amphtudes and timescales. Panels (B) through 
(H) of Figure [4] display large amplitude variable stars' 
light curves. Low amplitude variable stars are exhibited 
in Figure [5l The light curve featured in panels (C) and 
D) is one example of the fact that Kepler achieves a pre- 
cision sufficient to find transit-like features on the order 
of 200 ppm "by eye". 

5. CONCLUSIONS 

Kepler's grand voyage of discovery has begun and in 
the short while since Science Operations commenced we 
are already obtaining the photometric precision so in- 
dispensable to the task of discovering Earth-like plan- 
ets transiting solar-like stars. There are a number of 
instrumental systematics driven by changes in the tem- 
perature of the telescope and spacecraft that are both 
extrinsic to the photometer, such as the seasonal varia- 
tions due to the heliocentric orbit and quarterly rolls, and 
intrinsic to hardware on the spacecraft, such as the re- 
action wheel heater cycles. We are adjusting operational 



parameters in order to eliminate or reduce the ampli- 
tude of intrinsic systematics, or to move the timescale 
for these instrumental signatures outside of those rele- 
vant to transit detection, namely from 2 to 16 hours. 
These systematic errors represent nuisances that can be 
filtered out by data processing at the expense of some 
non-exoplanet astrophysics. The light curves obtained 
to date by Kepler are a delight to behold and promise 
unforeseen discoveries and a wealth of information about 
star's photometric behavior at a precision and over a con- 
tiguous and extended period of time not heretofore pos- 
sible. We look forward to establishing a tight constraint 
on the frequency of earth-like planets, a quantum leap 
in our understanding of the formation of planetary sys- 
tems, major strides from asteroseismic results and many 
serendipitous discoveries wholly unrelated to exoplanet 
science. 



Funding for this Discovery Mission is provided by 
NASA's Science Mission Directorate. We thank the 
thousands of people whose efforts made Kepler possible. 
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Fig. 1. — Scatter plot of the 0.5 hr-to-0.5 hr precision (A) and the 6 hr-to-6 hr precision (B) obtained in the Ql LC data set. The dashed 
and solid green curves are the upper and lower envelopes of the measurement uncertainties, respectively, propagated through the data 
processing steps used to construct the flux time series. Dwarf stars (log 5 > 4) appear as black points, while giants (log g < 4) appear as 
red points. There is a strong separation between the dwarfs and the giants in their photometric behavior at both timescales. Kepler is 
delivering near intrinsic measurement-limited noise performance across the full dynamic range of its target stars. 
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Fig. 2. — Time scries of changes in a) focus distance, b) mean row position and c) mean column position for the stars located on module 
2.1 near the outside corner of the FPA over a 10 day interval during Ql observations. 
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Fig. 3. — Raw and systematic error-corrected flux time series for one star on module 2.1 observed during Ql. The curves in panel A 
are the relative intensity of the raw flux (upper) and the systematic error-corrected flux time series (lower) for this star, which exhibits 
variability on timescales upwards of several days. Panel B displays the results of high pass filtering each of these time series to remove long 
term trends. The heater cycle pulses so prevalent in the raw light curve are suppressed in the systematic error-corrected light curve, and 
the long term variations are preserved since they are not correlated with instrumental systematic variables. 
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Fig. 4. — Relative flux variations for a variety of large amplitude variable stars and other flux time series exhibiting the broad dynamic 
range of Kepler's sensitivity and stellar variability. The left column displays the entire Ql light curve, while the right column shows close 
ups, except where otherwise noted. (A) Light curve for a star whose aperture is located on a column defect and "toggles" as its image 
drifts. (B) Close up on an RR-Lyrae star with an evolving envelope. (C and D) A flare star with flares that almost double the flux for 
this star on top of quasi-periodic variations of a few percent. (E and F) A star exhibiting high frequency oscillations on top of longer term 
quasi-periodic variations. (G and H) a highly eccentric eclipsing binary. 
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Fig. 5. — Relative flux variations for a variety of small amplitude variable stars and other flux time series exhibiting the broad dynamic 
range of Kepler's sensitivity and stellar variability. The left column displays the entire Ql light curve, while the right column shows 
close ups and zooms. (A and B) A star with high frequency oscillations on top of mmag variations. (C and D) Light curve for a star 
exhibiting 300 ppm transit-like features that can be seen easily by eye. (E and F) A star with variations well below 50 ppm, and coherent, 
small amplitude oscillations. (G and H) A quiet, main sequence star with a random fluctuations with an rms below 50 ppm on half-hour 
timescales. 



